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We consider the spectral statistics of the superposition of a random diagonal matrix and a GUE 
matrix. By means of two alternative superanalytic approaches, the coset method and the graded 
eigenvalue method, we derive the two-level correlation function X2(r) and the number variance 
E 2 (r). The graded eigenvalue approach leads to an expression for -^(r) which is valid for all values 
of the parameter A governing the strength of the GUE admixture on the unfolded scale. A new 
twofold integration representation is found which can be easily evaluated numerically. For A ^> 1 
qq ■ the Breit- Wigner width Ti measured in units of the mean level spacing D is much larger than unity. 

^\ \ In this limit, closed analytical expression for X2(r) and S 2 (r) can be derived by (i) evaluating 

^\ . the double integral perturbatively or (ii) an ab initio perturbative calculation employing the coset 

t-H ' method. The instructive comparison between both approaches reveals that random fluctuations of 

Ti manifest themselves in modifications of the spectral statistics. The energy scale which determines 
the deviation of the statistical properties from GUE behavior is given by VTl- This is rigorously 
shown and discussed in great detail. The Breit-Wigner Ti width itself governs the approach to 
qq , the Poisson limit for r — > oo. Our analytical findings are confirmed by numerical simulations of an 

■ ensemble of 500 x 500 matrices, which demonstrate the universal validity of our results after proper 

unfolding. 

t-H ! 

> . PACS numbers: 05.40.+j, 05.45.+b, 72.15.Rn 

oo ■ 

o\ : 

I. INTRODUCTION 

o ■ 

^\ , One of the archetypical problems in quantum mechanics consists of calculating (certain properties of) the eigenvalue 
""j—t ■ spectrum of a diagonal operator and a superimposed non-diagonal one. Little can be said in general about this 
problem. In our paper, we focus on the particular case where the matrix representations Hq and Hi of the above 
operators can be taken from the Poisson and the Gaussian Unitary EnsembleEl (GUE), respectively, 

H = Ho + aHi , (1.1) 

o 
o 



with a some strength parameter. The Poisson Ensemble is constructed from all those matrices whose eigenvalues 
are independent random numbers with identical, and largely arbitrary, distribution function. Due to the rotational 
invariance of the GUE, it suffices to consider only diagonal matrices Hq of the Poisson type. The random matrix model 



(1.1) should provide an adequate description for numerous transition phenomena from regular to chaotic fluctuation 
properties in atomic, nuclear, and condensed matter physics, as well as in quantum chaology (for a review see Ref. ||). 

Nuclear physics provides an important example. In heavy ion reactions, fast rotating compound nuclei are produced. 
The rotation is a collective motion of all nucleons which is often accompanied by single particle excitations. Thus 
the total Hamiltonian can be modeled as a sum of two contributions, a regular one, Hq, describing the collective 
motion, and a stochastic one, H±, describing the influence of the single particle excitations. This scenario was studied 
numerically in Ref. || for-.the cases that Ho is drawn from a Poisson or a harmonic oscillator ensemble. Using Efetov's 
supersymmetry methodQU, the qualitative behavior of the two-level correlations was also discussed analytically in 
Ref. |. 

As a further important ^xamplc we mention the problem of two interacting particles in a random potential in- 
troduced by ShepelyanskyQ. He predicted that two particles in a one-dimensional disordered chain can be extended 
on a scale L2 far exceeding the one-particle localization length L\. Subsequent wprkBlj quickly led to a definite 
confirmation and better understanding of this phenomenon. One possible approachO to this problem is to construct 
an effective Hamiltonian by diagonalizing the noninteracting part of the two-electron problem and expressing the mi- 
croscopic Hamiltonian in the basis of two-electron product states. The resulting representation consists of a diagonal 
contribution containing the eigenvalues of the noninteracting problem, and an off-diagonal contribution originating 
from the interaction operator. With the crucial assumption that both the diagonal and the off-diagonal matrix el- 
ements can be chosen to be random variables we arrive at the above random matrix model (for system sizes L of 
the order of L\). In the regime L > L\ the effective Hamiltonian has been studied in some detailllSril3. Imry has 
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shownB that the enhancement factor L%lL\ is given by the "Breit-Wigner width" or "two-particle Thouless energy" 
T2 measured in units of the two-particle level spacing D2, L2/L1 — T2. This raised, among other things, the ques tion 



how to identify 1-C2 in typical spectral observables like the number variance S 2 (r) of the random matrix model (1.1, 
considered herelij. This problem, and some surprises which we encountered while studying it, have been our original 
motivation in this project 



The random matrix model (Fl) has been considered by many authorstjOcJ'El. Recently, one of us succeedeoH 
in deriving exact integral representations. far the spectral fc-point functions of the Hamiltonian (F1L To this end, 
the graded eigenvalue method was usedlifrEj, which is a variant of Efetov's supersymmetry methodoO. lit, an effort 
to discuss certain approximations from these integral representations, two of the present authors derived 2 ^ a closed 
expression for the spectral two-point function X2 (r) by means of a special kind of saddle-point approximation. This 
led to the surprising observation that the energy scale at which -V 2 (r) and the number variance £ 2 (r) deviate from 
random matrix behavior is given by yTY and not by Ti. Here, Ti is the Breit-Wigner (or spreading) width induced 
by the perturbation aH\ and corresponds to the quantity T2 in the above example. Both T± andUr are measured in 
units of the mean level spacing D. This result was later confirmed in a perturbative calculationEil of the two-point 
correlator X2 (r) . The saddle-point approximation employed in Ref. ^ was limited to certain situations and does not 
include cases in which the energy separations r is much larger than the Breit-Wigner width Ti. 

The purpose of the present paper is threefold. First, we derive, avoiding the above-mentioned saddle-point expres- 
sion, a relatively simple exact expression for -^(r) for all r and all relative strengths of Ho and H\. For Ti ^> 1 this 
leads to very compact analytical formulas for both -^(r) and S 2 (r). The two scales y/Ti and Tx can be identified 
and interpreted. Second, we compare both variants of the supersymmetry formalism, the graded eigenvalue method 
and the coset method, by deriving our results independently for both methods. This sheds some light on the relation 
between these approaches. Furthermore our comparison elucidates the role played by statistical fluctuations of the 
Breit-Wigner width. Third, we confirm our result with the help of rather extensive numerical simulations of the 
spectral properties of the Hamiltonian (|1 . 1|) 



Our paper is organized accordingly. Following this introduction, Sec. O de als with the graded eigenvalue method 



the coset method is invoked 



and the improved treatment of the integral representations derived in Ref7|l8|. In Sec. Ill 
to essentially re-derive the results of Sec. O. We have tried to keep both sections reasonably self-contained. Readers 
who are completely unfamiliar with the supersym matr jy formalism should, however, consult the introductory literature 
for the cosetO'B and the graded eigenvalue methodEjIia. Section IV is devoted to our numerical simulations and Sec. |v| 
contains a summary and the discussion. 



II. GRADED EIGENVALUE METHOD 



In a first partqSec. II A, we introduce some basic terminology and briefly recapitulate the derivati on o f the integral 
representationslla of the fc-level correlation functions for the convenience of the reader. In Sec. II B the integral 



representation for the two-point function is further transformed to derive a relatively simple and, most importantly, 



tractable expression for the level-level correlator. Finally, Sec. II C is devoted to a perturbative evaluation of this 



expression in the limit where the Breit-Wigner width is much larger than the level spacing. 

A. Integral representation for spectral correlation functions 

The fc-level correlation functions for the Hamiltonian H in Eq. ([O]) characterize the spectral statistics completely 
and are defined by 

k 

R k (x u ...,x k ,a) = \ f P N (H)f[lmtr _ - - d[H] . (2.1) 
n * P =i x p 

The probability of finding k energy eigenvalues in infinitesimal intervals dxi around xi (i = 1, ...,k) is given by 
Rk(xi, . . . , Xk)dxi . . . dxk- For the case considered in this paper the probability distribution function P/v(i?), which 
depends explicitly on the dimension N of the matrices H, is given by the product 



P N (H)=pW(Ho)p£\H 1 ) (2.2) 



with 



2 



N 

n—l n>m 



oJV(W-l)/2 

4 1} (gi)= nN y 2 exp(-tr(ff 2 )) . (2.3) 



The function p^{H™ 1 ) is smooth but otherwise arbitrary. We introduce the modified fc-point correlators 
R k {x\, ... ,x k ), which are obtained by omitting the projection onto the imaginary part in Eq. ( |2.l| ). The original 
quantities can be reconstructed by appropriate linear combinations of the R k . To perform the ensemble average we 
write the modified correlators in terms of a supersymmetric normalized generating functional, 



1 d k 

Rk(xi, ...,x k ,a) = jTTTkTn 7TT Zk ( x + J ' a ) 

(2tt) k oJ\ . . . oJ k 



(2.4) 



.7=0 



where the energies and source variables form diagonal 2k x 2k matrices according to x = diag(xi, x±, . . . , x k , x k ) and 
J = diag(— Ji, Ji, . . . , — Jfc, Jfc), respectively. The averaged functional takes the form 

Z k (x + J,a) = J d[H ]Pjp(H ) J d[a]Q k {a, ^detg- 1 ^ + J - a) ® 1 N - l 2k ® H ] , (2.5) 

where a is a 2k x 2k Hermitean supermatrix and Q k (a,a) = 2 fc ( fc ~_L) cxp(— trger 2 /a 2 ) is a normalized graded proba- 
bility density. The subsequent steps can be summarized as followstJ. The matrix x + J is shifted from the graded 
determinant to the graded probability density and the supermatrix a is diagonalized according to a — u~ 1 su, where 
s = diag(sn, isi2, ■ ■ ■ , Sfci, isk2)- The volume element can be rewritten as d[a] — B%(s)d[s]dp;(u) with Bk{s) the Jaco- 
bian (called Berezinian in this case) of the transformation. The non-trivial integration over the unitary diagonalizing 
supergroup with its Haar measure dfjb{u) is the central step in the graded eigenvalue method and can be performed 
with the supersymmetric extension of the Harish-Chandra Itzykson Zuber integral. Collecting everything we arrive 
at 

Z k (x + J.a) — 1 ~ rj(x + J) + 1 / G k {s-x- J ia )zf\ S )B k {s)d[s\ , 

±>k{x + J) J 

Zf{x + J) = J rf[H ]P^ 0) (i? )detg- 1 [( a; ± + J) <g> ljv — hk ® H ] , (2.6) 
where the kernel resulting from the group integration is Gaussian and given by 

G k (s -x,a)= 1 exp ( — \ trg(s - xf ] . (2.7) 

(7TQ^) K \ a z / 



The distribution 1 — rj(x + J) in Eq. (2.6) ensures the normalization Z k {x,a) = 1 at J = 0. We mention in passing 



that the generating function Zfc(izi + J, a) satisfies an exact diffusion equation in the curved space of the eigenvalues of 

Hermitean supermatrices. Here, t = a 2 /2 is the diffusion time apd-the generating function zj^\x + J) serves as the 
initial condition. This diffusion is the supersymmetric analogueliSE 2 ] for the diffusion of the probability distribution 
function in the space of ordinary matrices, which is equivalent to Dyson's Brownian. Performing the source term 



derivatives in Eq. (2.4) we find 



R k ( Xl ,...,x k ,a) = ^- f G k (s-x,a)$Z ( k 0) {s)B k ( S )d[s} , (2.8) 



where the symbol S3 indicates the above-mentioned proper linear combination of terms. In order to arrive at truly 
universal quantities we have to measure all energies in units of the mean level spacing D. This leads to the definitions 
i P = Xp/D, X = a/D, Xk(£i, ■ ■ ■ ,ffc) A) = hin^^oo D k R k {xi, . . . , x k , a), and z£ 0) = limAr^oo zf ] (Ds). The level 
correlation functions X k defined in this way are translation invariant over the spectrum. They can be expressed as 



X k (£ 1 ,...,t k ,\) = L-9- [ G fe ( S -^A)3z[ 0) ( S )i? fc ( S M S ] 



(2.9) 

after the redefinition s — > s/D. We notice the very similar structure of the integral representations ( |2.8| ) and ( p.9| ). 
This is so because the above mentioned diffusion in superspace is, apart from the initial conditions, the same on all 
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scales, in contrast to the diffusion in ordinary space which is equivalent to Dyson's Brownian. The two-point function 
X 2 , to which we res tric t our attention henceforth, depends on r — £1 — £2 only and we can perfawn two of the four 
integrations in Eq. (2.9) due to translational invariance. This leads to the integral representationtS 



X 2 (r,A) = - 



7T 3 A 2 



OO OO 

/ / ° XP (~2^ ( ** + ^ } ) ^ k %^% W+^y> 9f 4° ) (*i,«a)'fti<fta 



(2.10) 



—00 —00 



with ti = sn — S21 and t 2 — s 12 — s 2 2- The initial condition is still arbitrary. For the case of Poisson regularity 
we have 



$szi°\ti,t 2 ) = -Re ( exp ( -ztt n ^ - 1 



1 



2it; 



(2-11) 



Equation (2.10) is exact but, unfort unate ly, difficult to evaluate as it stands. In the following subsection, we derive a 
more convenient formulation of Eq. ( 2.10] ) . 



B. Simplification of the two— point level correlation function X%(r, A) 



As explained in the introduction, a previous to extract physical information from Eq. ( |2.10| ) led to the 

discovery that the energy scale at which X 2 (r, A) deviates from random matrix behavior is linear in A. For very 
general reasons, however, the Breit-Wigner width is always quadratic in the strength of the perturbing matrix 
elements, hence T± cx A 2 . This led to the conclusion that y/Ti is the important energy scale. The approximation used 
to derive this result was, however, not valid for r 3> IV Here, we avoid the saddle-point approximation proposed in 



Rcf. 



and proceed as follows. We introduce the new variables Xj 

1 r 

7rA^ 



tj/(Tr\ ) and the important abbreviations 



^ 2 A 2 ' 



(2.12) 



It is instructive to note that ir\ 2 is actually the Breit-Wigner width measured iu units of the level spacing, 7rA 2 = IV 
This can be deduced from the local density of states for the model studied herd 2 !! and will become more obvious when 
we discuss the coset method in the next section. With the above notation we find 



X 2 (r,X) = X 2 (k,c) = X 2 (k,c) + X 2 (-k,c) + n 2 



where 



(2.13) 



X 2 (K,C) = 

7T 



dx±dx 2 - 



xix 2 



—00 —00 



C(xi,x 2 ) 



1 



X\ + i 



(xi 



x\+xlY 



exp(-£(xi, x 2 )) 



— (X 2 + IK —J 



Xi+l' 



X 2 {k,c) 
£-{x x ,y 2 ) 



dxidy; 



(xi + i)(y 2 + in) 



1 



-OO —OO 



2c x\ + i 



(xi + i) 2 + y 2 
[(xi + i + n) 2 + (y 2 + in) 



exp(-£(xx,y 2 )) 



(2.14) 



To reduce the order of the pre-exponential singularity we perform an integration by parts using 2x 2 /(x1 + x 2 ) 2 = 
—d X2 (x 2 + x 2 )^ 1 . With the additional transformation x 2 = y 2 x\j(x\ + i) Eq. ( p. 14 ) takes the form 



(2.15) 



Now, our strategy is to eliminate the imaginary contributions to the squares in C(xi, y 2 ). To this end we perform two 
shifts, namely y\ = x\ + % and u 2 = y 2 + in. This amounts to moving the integration contour from the real axis to 
the lines Im(xi) = —1 and lm{y 2 ) = —k, respectively, see Fig. |l|. The integral (2.15) remains unchanged under these 
shifts unless the integrand exhibits singularities in the region between the old and the new integration contours. With 
the help of Fig. [l]the reader can easily convince himself that such singularities do indeed exist in our case. They give 
rise to two different residuum contributions R\ and R 2 . We therefore arrive at 
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X 2 {k,c) 



dyidu^ 



yiu-2 



])\ + («2 - in)' 



■ exp(— C(yx - i, u 2 - in)) - Rx - R 2 



Ri 



i?2 



dy 2 {V2 + J dy 2 (y 2 + in) exp ( 2 «^(1 - 2/2)) 





K + \k\ 



dy 



yi - k) exp ~{yi - i) 



dyi 



y x + k exp -(yx - i 



K — \k\ 



(2.16) 



We notice th e app earance of sign.(/«) = k/\k\ which is needed to distinguish between the two cases k < and 
K > in Eq. ( 2.13 ). The residuum contributions can be easily calculated analytically and with the fin al (in nocuous) 
transformation ux — y\ + k we obtain the following result for the two-point level correlation function (2.13), 



X 2 {r,X) = 1 + 



2(nr) 



2 

00 00 



exp 



J A 2 



duxdu 2 



cos(27rr) — 1 
(ux — k)u 2 



(ux — n) 2 + (u 2 - in)' 



■ exp 



1 

27 



ii\ — n 



(uf+ui) +{k 



(2.17) 



In the limit A — * 00, the first two terms in Eq. (2.17) reduce to the two-point level correlation function for the GUE, 

2 



X 



GUE 



(r) = l 



sm 7rr 



(2.18) 



The additional terms are corrections for finite A. We denote the double integral in Eq. (2.17) by P(c,k). It can be 
evaluated numerically for a wide range of values of c and n without too much effort. For our present purposes, in 
particular for the comparison with the coset method in Sec. Ill, we prefer to derive some analytical results in the 
limit c«l. 



C. Perturbation theory for X2(r,X) 



In the following we assume that c — > 0, but k — const. This means that the Breit-Wigner width is much larger than 
the level spacing, I"i 3> 1, and that the energy difference r scales like Ti. Under these circumstances P(c,k) can be 
calculated perturbatively, resulting in a power series in c. With the polar coordinates ux — p sin ip and u 2 = pcosip 
and after expanding both the exponent and the pre-exponential term in Eq. (2.17) in powers of p we obtain the 
following expression for the double integral, 



P(c, k) — — dip dp pe ltp cos ip ( 1 — — sin <p ) ~S~" ( ^— ) exp 



m=0 



x exp 



1 K + i 2 

2 CK 



2ck 



E 

n=l 



psmip 



(2.19) 



The structure of the angular, i.e. the ip, integration indicates that only terms containing odd powers of p lead to 
non- vanishing contributions to P(c,n). From Appendix [a| it is clear that P(c,k) comprises terms of order c 1 , c 2 , 
c 3 , and so forth. For later reference, we will calculate the c 1 and c 2 contribution explicitly. To first order in c the 
contribution of the double integral P(c, k) to the two-point correlator X 2 (r, A) is (cf. Appendix [A|) 



P^(c,n)+P {1 \c,- K ) = 



ck ck 1 Tx f r 

2{k + i) + 2{k - i) = nr 2 +Tj \ Y~x 



(2.20) 



where we have used that Tx = n\ . This result can be combined with the term 1/(tt X ) — l/(7rr 1 ) = c in Eq. (2.17) 



to give the full first-order contribution X^\r, X) to X 2 (r, X) 
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X«(r,A) = - 



i ri 



7T r 2 + r( 



(2.21) 



The se cond -order contribution X\ (r, A) is composed of four different combinations of terms from the perturbation 
series (2.1£), see Appendix [A|. Combined they can be written as 



X. 



.2) 



1 



2/ 



1 



r?- 



2tt 2 (r 2 + T 2 ) 2 



(« + i) E 
c 2 
4 



(K + i) 4 

2i 



+ (k «-> — k) 
3 



(«• 



+ (k «-> — k) 



(2.22) 



(2) 

In the last line of the above equation we have distinguished between two contributions to X\ (r, A) , a first piece 
originating from the term with the quadratic denominator, and "extra" contributions arising from the terms with 
higher-order denominators. The significance of this distinction will become clear in the following section, where we 
essentially re-derive the present results from the point of view of the coset method. 



III. COSET METHOD 

The supersymmetric "coset" methodBi has been widely employed in the last_ten years or so to solve problems 
of random matrix theorya. In the present context, it was already used in Ref.0 for a qualitative discussion. To 
make contact with the well-established methods used in the literature over the years we will adopt normalization 
conventions in this section which differ slightly from those used in Sec. ||. 



A. Basic definitions 



We write the full Hamiltonian H in Eq. (1.1) as 



H = H + H 1 



(3.1) 



thereby effectively absorbing the strength parameter a into the definition of Hi . The distribution (77) of the non- 
vanishing elements Hq of Ho defines an important energy scale, the bandwidth B, because we have p( 0) (Htf) oc B^ 1 
for reas ons of normalization. The probability distribution function pffi (Hi) adopted here is a slightly modified version 
of Eq. (2J3) and reads 



oil) 

N 



(Hi) = 2^-D/ 2 f 



N 
2tty 



N 2 /2 



exp 



N 
2o~ 



(3.2) 



where the symbol (. . .) denotes averaging over the _ffi-ensemble (i.e. the GUE). The parameter 7 introduced in this 
way is related to the strength parameter A employed in Sec. [n] through 



A 



N D 



(3.3) 



It is our purpos e in this section to calculate the two-point correlator X2 (and later the number variance S 2 ) for 
the Hamiltonian (3.1) perturbatively in a certain suitable range of parameters. To define this range we note that 
D = B/N is the level spacing of Hq and j/VN the typical strength of the perturbing matrix elements. Hence we 
find for the (dimensionlcss) induced spreading or Breit-Wigner width Ti 



TiD oc 



ND 



B 



(3.4) 



Our calculation in this section is valid under the two conditions 



1_ 

B 



TiD 
B 



Ti oc -!— > 1 
1 BD 



(but finite) . 



(3.5) 



G 



The first of the conditions (3.5) means that the bandwidth B is infinitely larger than the Breit-Wigner width T\D. 
This ensures that neither B nor the level spacing D are appreciably changed by the perturbation H\. The second 
condition ensures that we deal with overlapping resonances, i.e. the original states of H are thoroughly mixed to form 
the new eigenstates of the combined Hamiltonian H. Our perturbative calculation relies on the small parameter 1/Ti. 
Again, as in the previous section, the parameter n = r/T± is held fixed so that the dimcnsionlcss energy difference r 
scales like Ti, 

With the Green function G± — (x ± — H)^ 1 the density of states is given by 

P (*) = ^(trG+-trG_) (3.6) 

and hence the averaged density-density correlation function takes the form 
p 2 (x + u/2, x - w/2) = [(p (x + lu/2) p(x~ w/2) )] „ 

= -^Re[(trG+(a; + w/2)(trG_ (x — u>/2) — tr G+ (x — w/2)) )] Q . (3.7) 

Here, we have introduced the notation [. . ,]q for the i?o _ average. The angular brackets (. . .) denote, as above, averaging 
over Hi. In the sequel we will typically perform the -Hi-average in a rather early step of the calculation while the 
-ffo-average is performed in the final stages. We note that 



/M|,--hu.-/2..,--u.-/2) = ±X 2 (r) + ^8(r) (3.8) 



defines the relation of the density-density correlator to the spectral two-point function (2.9) considered earlier. Here 
r = lo/ D is the dimensionless energy difference as in Sec. 0. 



B. One— point function and density of states 



To illustrate our procedure and to define relevant energy scales we start with the average one-point function. With 
the supersymmetric generating functional 



Z±{q±) = detg -1 ^ - H] ® 1 2 + ljv 8> q±L g /2) 
= exp (-trgln (ljv <g> 1 2 + G± ® q±L g /2)) 



we can write the Green function as 



trG±(i) 



d 



d, 



Z±(q±) 



q± 



(3.9) 



(3.10) 



q=0 



where L g = diag(— 1,1) is a two-dimensional (super-) matrix. With standard techniquesoO the £/i-average of the 
generating functional can be easily calculated to give 



(Z+(jj) = J d[a] exp f-ytrg(a 2 ) - trgln([a; + - H ] ® 1 2 + In ® [ja + J}) 



(3.11) 



We have introduced the source matrix J = qL g /2, and a is a 2 x 2 supermatrix as in Eq. (2.5). In contradistinction 
to the previous section the further evaluation of Eq. (3.11) relies on the saddle-point approximation. With the 
saddle-point equation 



N ^ 



1 



, x+ - Hi 3 + 7 cr 
and with the ansatz 717 = (To + iT±)D/2 = TD/2 we obtain 

2 7 2 ^ 1 



(3.12) 



r = — 



ND 



2f_ 
D 



E 



x+ - Hi 3 + DT/2 



p(°)(r/) 
x+ + DT/2 



(3.13) 
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Equation (3.13) can be solved in the limit |T| <C B and we arrive at the approximate expressions 

r 



.?p/d/ (0) W 



D 

As immediate consequences we find 



D j x — j] 
r 1 «27r-^„W(0)=27r^ 



(trG+) = tr- 



1 



(Breit-Wigner width) . 



NTD 



(3.14) 



x+ -Hq + DT/2 



2 7 2 



(p( 3 ;)>=-ilm<trG + > 1 = ^=^)(x) 



(3.15) 



In particular, the last equation means that the Hq average of the mean level density can be trivially performed, giving 
[(p(x))} = Np^(x) 

We see from Eq. ( |3.7| ) that we need the average of a G + G_ and of a G + G + term. It is well known that in the large 
N limit the H\ -average of the product of two Green functions with infinitesimal increments of equal sign factorizes, 



(trG + (x + uj/2)trG + (x-uj/2)) = (trG+)(trG+) = tr 



-tr- 



x+ + u>/2- Hq+DT/2 x+ -u/2- Hq + DT/2 



(3.16) 



The same turns out to be true after averaging over Hq. This average is easily performed by simply integrating over 
the independent random entries of the diagonal matrix Hq, 



[(trG+trG+>] «iV(iV-l) / dr, 
+N / dr] 



u/2-T] + DT/2 



dr) 



p (0) (v) 



x+ - uj/2-r] + DT/2 



(x+ +u/2-r] + DT/2) (x+ -u)/2-r) + DT/2) 
^^[(trG + )] [(trG + )] « [(trG + )] [(trG + )] . 



(3.17) 



The term oc N in Eq. ( |3.17 ) is seen to vanish upon closing the contour in the lower half plane. We conclude that in 
order to obtain the density -density correlator (3.7) we can simply replace (trG±) = +vk/D. 



C. Two— point function and nonlinear a model 



Our next goal is to calculate the remaining quantity (trG+trG-), which leads to the familiar no nlinea r a model, 
and its -ffo _ average. First, we generalize the definition of the supersymmetric generating functional ( |3.11 ), 



(tr G + (x + c/2)trG_(x - w/2)> - — — (Z(q + , q_)) 

Oq + Oq_ 



q±=0 



d d 

d q+ d q _ 



d[Q]exp(-£(Q)) 



q±=o 



C(Q) = J2 tr § ln ( N - H o'} U + | A + ^TQ + J 



(3.18) 



The quantity Q is a 4 x 4 supermatrix parameterizing the saddle-point manifoldala to which the integration has been 
restricted by the saddle-point approximation. Furthermore, J = diag(<7+L g , q-L g )/2 = q+ P+ + q~P- (this defines 
the projectors P + and P_), and A = diag(l2, — 12)- To perform the derivatives in Eq. ( [3.18 ) it is useful to expand C 
up to second order in the source matrix J, 



C(Q) = Cq{Q) + d{Q) + C 2 {Q) 
The three contributions can be written as 



(3.19) 
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Co = trgln[(z - Hi j )l A + (wA + r D + iT 1 DQ)/2] 

3 
3 

£■2 = -l^2^g[93 J 9jJ] , 



(3.20) 



where we have introduced the abbr eviation gj — [x + (w/2)A — Hq 3 + (TqD + iTiDQ) /2] l . After performing the 
source term derivatives in Eq. (B.1S) we arrive at a central equation of this section, 



(trG + (x + w/2)trG_(a;-w/2)> = J d[Q](5 x (Q) + S 2 (Q)) exp(-£o(<9)) 

S 1 (Q)=^trgfeP + ]trg[ fffe P_] 



S 2 (Q) =J2^[9jP+9 3 P-} 



(3.21) 



This is the particular form of the zero-dimensional nonlinear a model describing the crossover between Poisson and 
GUE statistics. 



D. Evaluation of the nonlinear a model 



As already mentioned above, the Q-integration in Eq. ( 3^21 ) is restricted to the saddle-point manifold familiar from 
numerous previous applications of the nonlinear a modelH It belongs to the peculiar features of superanalysis that 
integrals of the type (3.21), which are derived by a supersymmetric change of variables, contain an extra "boundary" 



contribution (sometimes referred to as the Efetov-Wegner term). This boundary contribution is generically given by 
the value of the integrand at Q = A. In our present case we have £q{A) = 0, 52(A) = 0, and 



1 



j^x- Hi 3 + DTa/2 + (uj + iT x D)/2 x - H* k + DT /2 - (w + iT 1 D)/2 
The i?o~average is performed as in Eq. ( |3.17| ) , 



(3.22) 



Si (A) = N(N — 1) / dn 
+N / d-q 



x-n + DTq/2 + (lu + iTiD) /2 



d.rj 



x-n + DT /2 - (u + ?TiD) /2 



(x - r? + DTq/2 +(uj + iT x D)/2) (x - n + DT /2 - (w + iT x D)/2) 

If, 1 



[(tr G+ (a; + u/2))] Q [(tr G_ (a; - w/2))] Q + ^ 



z 2 - {uj + iT^Y/A 



= [(tr G+ (rr + w/2))] Q [(tr G_ (z - w/2))] Q + 



2tt 



1 



Z? 2 Ti — ir 



(3.23) 



We recall that r = lo/D. The first term in Eq. (3.23 together with Eq. (|3.17| ) for ms th e disconnected contribution to 
the two-point correlation function A^r). Upon inserting the results ( ft.l4| ) and ( 3.1 5| ) one can e asily show that the 
sum of disconnected contributions to A^r) reduces to unity (for To « 0). The second term in Eq. ( 3.23| ) is part of the 
connected contributions. In the following we calculate the remaining connected terms by treating the saddle-point 
integration perturbat ively. 
We start from Eq. 



(3.21 



The general strategy will be to express both the exponent and the pre-exponential terms 
in Eq. ( [3.21 ) in terms of the independent variables of the saddle-point manifold. It turns out that we have to consider 
only quadratic terms in the exponent so that the integration becomes Gaussian and therefore trivial. Concerning the 
-ffo - & v erage it is important to note that those terms involving correlations between the matrix elements of Hq in the 
exponent and in the pre-exponential terms can be neglected. In fact, it turns out that the pre-exponential terms and 
the exponent can be independently averaged over Hq. This amounts to a tremendous simplification of our calculation. 
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A more thorough discussion of these issues as well as a number of technical steps omitted here for clarity can be found 
in Appendix |b|. 

We express the deviation of Q from the diagonal value A in terms of the quantity AQ — {A, Q}/2 — 1, where {., .} 
denotes the anticommutator. Then the i?o~average of Cq(Q) can be expressed as (see Appendix (|B|)) 



[ £ o(Q)] --^ ! f £ ^ trg(Ag) - 

For our present purposes, a suitable parameterization of the saddle-point manifold is given by 

Q = T _1 Ar = AT 2 



(3.24) 



T= \/l + R 2 + R 
t 
t 



R 



(3.25) 



where t and t are 2x2 supermatrices representing the unrestricted "free" variables of the saddle-point manifold. It 
follows for Co that 



[MQ)} 



-2wi J" \ trg(tt) 

I i — IT 



(3.26) 



Likewise we have to express the -ffo _ averages of S\ and 1S2 in terms of t and t. Again, the detailed derivation has been 
deferred to Appendix H, 



[Si 



rj 



D 2 (Ti 



r trg[i7P + ]trg[ttP_] 



(3.27) 



D 2 (T x - ir f 

Using Eqs. j3~2l| ), (j|J|), and ( ^2?| ) the remaining Gaussian integrations can be easily performed. To this end it is 
useful to employ certain Wick-type contraction rules as explained in Ref. |2^. Along these lines we obtain the following 
"diffusion" contribution from a perturbative treatment of the saddle-point integral, 

r? 



[<trG + trG_)] 



diff 



1 



2il\ 



1 



1 



(Dr) 2 (Ti - ir) 2 D 2 r{T 1 - ir) 2 (Dr) 2 D 2 (T 1 - ir) 2 



(3.28) 



The individual contributions of Si and 1S2 to this result correspond to the first and the second term on the r.h.s. of 
the first line, respectively. 

In summary, the t otal (Ho and Hi) average of trG+trG- consists of three terms, the disconnected part, the 
boundary term ( 3.23| ), and the diffusion contribution ( p. 28] ), 



[(trG + trG_)] = G disc + G, 



bound 



Gdiff ■ 



By definition, see Eqs. (p.7|) and (p.8[), we have 



D 



X 2 (r) = ^Re(G boU nd + C diS ) . 



(3.29) 



(3.30) 



If we insert our results for Gbound and Gdiff and generalize the latter to arbitrary (GUE, GOE, GSE) symmetry by 
multiplying it with 2/(3 (j3 = 1,2,4) we arrive at 



X 2 (r) = l-i 7 4r 5 

2 ' it r? -L- - 2 



(2), 



4 2) (r) 



1 T\ 



(3.31) 



tt 2 /3 (T| +r 2 ) 2 ' 

The first term in this result is the perturbative (divergent) expression for the two-point co rrela tion function of the 
Gaussian ensembles. It can be replaced by the full non-perturbative result, i.e. by Eq. ( 2.18 ) in the case of the 
GUE. Comparison of the first and second order contributions to A 2 (r), X^\r) and X% (r), obtained so far with the 
corresponding results of the graded eigenvalue method in Eqs. (2.21) and (2.22), reveals missing terms in the present 
calculation. We conclude this section with a discussion of the origin of this discrepancy. 
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E. Fluctuating Breit— Wigner width 



Let us recall the saddle-point equation ( 3.13 ). From the first line of this equation it is clear that T depends on 
the particular realization of Hq. But so far we have only considered one fixed (mean) value of T and have neglected 
random fluctuations of T within the ensemble. We will see below that these random fluctuations of the Breit- Wigner 
width lead to the missing terms in our expression for Xv(ui). To discuss the effect of fluctuations we choose the 
distribution of the (diagonal) matrix elements of Hq to be 



1 

2W 



e(w-\ v \) 



(3.32) 



for definitcncss. This means that the bandwidth B is given by B = 2W. Then, according to Eq. (3.14), 

r -o 

27TT 2 

Ti = Nh {D = 2W/N) ■ 

We note that if we insert the rel ation ( |3.3| ) between 7 and A in Eq. ( 3.33| ) we obtain Ti = tt\ 2 as claimed below 
Eq. ( 2.12 ). In the approximation ( 3.33 ) T is given by 



(3.33) 



r<°> = . 

With this notation the exact saddle-point equation fl3.13| ) reads 

IYD 



r = 



£ 



x - Hl° + DT/2 



(3.34) 



(3.35) 



If we insert our lowest (zeroth) order approximation (3.34) on the r.h.s. of Eq. (3.35) we obtain the next (first) order 
in an iterative refinement procedure, 



r 



(i) _ 



TiD 



£ 



1 



jy« +jDr (o)/2 ' 



so that the random fluctuations of T c are given by 

d T = r (1) -r (0) =-^— 



E 



x - Htf + iDY 1 /2 



(3.36) 



(3.37) 



These fluctu ations lead to additional terms in the boundary contribution because their inclusion modifies the r.h.s. 
of Eq. d3~22] ). With the definitions 



G, 



x + u/2- Htf +L>r(°)/2 
1 

x- lu/2- H[ j +DT(°>/2 



(3.38) 



we have 



x + lu/2- Hi 3 + D(T(°) + ST)/ 2 



x-uj/2- Hi 3 + D(r(°)* + ST*)/2 



G . - E-6TG* 



— D — 2 

Gj ST*G, 

o 2 3 



(3.39) 



Hence we obtain instead of Eq. ( 3.22 ) 

S, (A) = ^ (g s - f STG^ ^ (G k f ff*G 2 fe ) 

st = -^J2g 1 -*t 1 

I 

7T ' J 



(3.40) 
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Averaging Si (A) over Hq is tedious but in princi ple st raightforward (see Appendix The result of this calculation 
is that there are additional contributions to Eq. ( 3.23 ), 



/-textra 
^bound 



2r a 



3r 2 



D 2 {T l -irY D 2 {T l 



leading to the extra terms 



X. 



(2) 



c 

T 



2 1 



(K + i) 3 (n + i)' 



+ (k «-> — k) 



(3.41) 



(3.42) 



in the second order contribution to X 2 (u;). These are precisely the missing terms in Eq. (3.31). We conclude that 
random fluctuations of the Breit-Wigner width manifest themselves in the two-point level correlation function. These 
extra contributions are automatically included in the exact graded eigenvalue method. 



IV. NUMBER VARIANCE AND NUMERICAL SIMULATIONS 

The number variance £ 2 (r) for a random matrix ensemble is defined as the variance of the number of energy levels 
within a dimensionless energy interval of fixed length r. It can be easily calculated from the spectral two-point 
correlation function X 2 (r, A) , 



£ 2 (r) = r - 2 / (r - r')(l - X 2 (r', X))dr' 


r 

= -r(r - 1) + 2 J(r- r')X 2 (r', X)dr' 




(4.1) 



This spect ral observable is particularly suitable to identify the energy scale at which the spectral statistics of the 
ensemble (1.1) deviates from either the GUE or Poisson behavior. For the Poisson ensemble, £ 2 (r) is simply given 
by the straight line, E| oi (r) = r, while for the GUE t he n umber variance is of roughly logarithmic shape. For 
intermediate cases, i.e. for finite non- vanishing a in Eq. (1.1), one expects £ 2 (r) to exhibit GUE behavior for small 
energy separations r. For lar ger e nergy intervals it should the cross over to a linear Poisson-like form. In Sec. IV A 
we determine, based on Eq. ( ff.l| ) a nd th e correlator X 2 (r,X) calculated in the previous sections, the energy scales 
relevant for this crossover. Section IV B is then devoted to the question whether our theoretical expectations are 
borne out by direct numerical simulations. 



A. Crossover Energy Scales 



The principal result of the analysis performed in Sees. || and 111 for the case Ti S> 1 is that the two-point level 
correlator X 2 (r, A) for the random matrix model flL.l| ) can be written as 



X 2 (r, A) = X 2 GUE (r) + X 2 NA (r, A) + X^ (r, A) + X^ (r, A) 



(4.2) 



* 2 GUE (r) = 1 - 



X 2 NA (r,A) 



xP(r,X) 



sin 2 (nr) 
(irr) 2 



1 



2tt 2 r 2 

1 r x 



= 1 - 



1 



2tt 2 r 2 
.2 



exp — 2-7T- 



7T T 2 + r 2 

1 T\-r 2 



1 



(1 — cos(2-7rr)) 
cos(27rr) 



27T 2 (r 2 +r 2 ) 2 2tt 2 V (r 2 + r 2 ) 3 



1 ^-2r 2 (r 2 -3r 2 ) 3r 2 (rf-6r 2 r 2 + r 4 ) 
f (r 2 + r 2 ) 4 
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with Ti = 7rA 2 = n(a/D) 2 . The first contribution Xp UE (r) is just the two-point level correlator for the GUE case. 
The contribution X^ A (r, A) describes nonanalytic corrections for T^ 1 <C 1 (with r of order Ti) that cannot be obtained 
by perturbation theory. They are found in the framework of the graded eigenvalue method or a more sophisticated 
nonperturbative evaluation of the coset method. This contribution tends to eliminate the oscillating part of X 2 UE (r) 
for r ^> \Zr~i. The terms X^\r, A) and X^(r,X) are just the first or second order corrections in T^ 1 . In (E73), 
we have focused on the unitary case {(3 = 2) for definitcness. Our qualitative conclusions, however, apply likewise 
to the orthogonal and the symplectic symmetry class (/3 = 1, 4). The nonanalytic correction can be neglected when 



performing the r'-integral in (4.1) to determine the number variance. As a consequence of the form of Eq. ( [1.2] ) we 
obtain 

j: 2 (r,X) = E 2 GVB (r) + E 2 {1) (r,X) + E 2 {2) (r,X) . (4 . 3) 



A short calculation using Eqs. (4.1) and (4.2) reveals that 



2r / r \ ( T 1 1 \ / ( r \ \ 1 r T i 



A*(r, A) = £ 2 1} (r, A) + £ 2 2) (r, A) = -arctan ( - ) - ( -± + ^ ) In ( 1 + [ K ) ) + ^T^TtW ' ( 44 > 



T X J \tt 2k 2 J \ V r i/ / 27r 2 (r 2 +r 2 ) 2 



It is consistent with our assumption Ti ^> 1 to neglect the second term l/2n in the prefactor of the logarithmic 



contribution and the last term in Eq. (4.4). To identify the energy scales of interest we first consider the limit r<Ti, 
Then 

Obviously, the scale set for deviations from the GUE behavior is \fT\ in agreement with previous findingsll. 



Let us now ask a different questio n, n amely, what is the limit of Eq. (4.4) for r T± and how is it approached? 



We expand the arcustangent in Eq. (4.4) for large argument and arrive at 




AS 2 (r,A)«r--r 1 -^ln 1+ — . (4.6) 



Apar t fr om a constant and logarithmic corrections in r (in this context the GUE result T,Q VE (r) to be added to 



Eq. ( [4.6| ) is also nothing but a logarithmic correction) we recover the Poisson limit Ep oi (r) = r. The en ergy scale 




governing the asymptotic approach to the Poisson limit is given by T\. We can see this by dividing Eq. (4.4) by r. 
Within our usual approximations we find 

J-AE 2 (r, A) = ^arctan ^ ) - ^ In | 1 + ( — ) ) . (4.7) 

This function approaches 1 for r ^> Ti and depends only on the ratio r/Y\. Hence the relevant energy scale is Ti as 
claimed. 

It is quite instructive to compare our results with the case of a disordered metal in the diffusive limit. The two- 
point level correlator for this case was calculated by Al'tshuler and ShklovskiiEa in the framework of diagrammatic 
perturbation theory, 

X 2 = 1 + -^r Re V — ^ „.„ (4.8) 

\DJ /3tt 2 {~iuj + x> q 2 ) 2 

where (3=1, 2, 4 labels the symmetry class and T>o is the diffusion constant (note that D is the mean level spacing). 
This expression is valid for lo 3> D. The sum extends over the diffuson modes characterized by wave vectors q. The 
mode q = corresponds to random matrix theory while the modes q ^ give the deviations from random matrix 
behavior. The latter can be described by a function f(x) in the following way 

M£) =1 -^ + f '(£)< /M=«oW + .... (4.9) 

Here E c ~ Vq/L 2 is the diffusive Thouless energy for a sample of size L. The coefficients clq, a 2 , ■ ■ ■ and therefore the 
function f(x) depend only on the space dimension d and on the particular geometry of the sample but not on any 
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other physical parameter. These coefficients are given by discrete sums over "integer" wave vectors which converge 
nicely for d < 3. The correction of the nonzero modes to the number variance is immediately estimated as 

A ^~ww + --- (4 - 10) 



with r = uj/D. This corrections is of order 1 just for r ~ E c /D and not y/E c /D. Here the crossover is indeed governed 
by only one parameter and there is no transition regime like the interval between %/rT and Ti as found above. 



B. Comparison with Numerical Spectra 



For a verification of our analytical calculations we numerically generate a matrix ensemble of the type (1.1). The 
distribution function (r)) of the independent elements of the diagonal matrix Ho is chosen to be a semicircle with 
radius Rq = 30 (in arbitrary energy units). The matrix Hq is perturbed by a GUE ensemble Hi with (\H l ^\ 2 ) — a 2 /2 



as in Eq. ( |1 . 1[) . We consider matrices of dimension N = 500 and for each set of values a, Ro we take 500 realizations 
into account. To ensure that I\ > D in our simulation we take the parameter a from the interval [0.1,3.0]. Larger 
values of a are no longer useful because it is difficult to distinguish the ensuing spectral statistics from pure GUE 
statistics. 

The spectra obtained by diagonalization of H = Hq + Hi are unfolded numerically. Then the number variance 
E;? im (r) is calculated for a number of diffe rent v alue s of a. Given this set of numerical functions we fit (by hand) 



the theoretical formula given by Eqs. (4.3) and (4.4) to the data. This allows to determine the unfolded crossover 
parameter A (via the relation I\ = tt\ 2 ). Furthermore, we also extract the "critical" energy scale r c at which £j? im (r) 
deviates from GUE behavior from the data. We define r c by the condition 

^ m ^l ~ ^ue(^) = , , , 4 , | , 



which means that at the point r c the actual (simulated or fitted) number variance exceeds SQ UE (r c ) by 20 percent. 
Of course, this condition is somewhat arbitrary. However, in our opinion it correctly reflects what we usually mean 
by a "significant deviation from GUE behavior" . By means of the procedure just described a given value of a is 
associated with particular values of A and r c . The following figures illustrate the procedure and our results. 

In Fig. ^|we show the numerically calculated number variance £g im (r) and the corresponding fit for the two values 
a = 0.5 and a = 1.2. In the latter case, there is very satisfactory agreement up to very large energy separations 
of r = 140. In the former case, we find agreement of this quality up to r = 30. Beyond this value, we find some 
disagreement which is likely to be due to the limits of our approximation which relies on large A values (i.e. a 3> D). 
Fig. H shows the fitted crossover parameter A and the theoretical value A t h = a/D. For small values a < 0.5 
both curves agree well. The deviations for a > 0.5 are due to the fact that the Breit-Wigner width I\-D already 
approaches the effective bandwidth of the unperturbed spectrum. In fact the radius Rq of the semicircle density equals 
TjD for a c » 0.95. In this regime both theoretical approaches arc strictly speaking not valid due to the condition 



(3.5). However, it seems that the theoretical expression for the number variance remains even here correct, if the 
parameter A is correspondingly rescaled ( "fitted" ) . Fig. [| shows the dependence of the critical energy separation r c 
on a. Two observations are important: First, the dependence is by no means linear. Second, the data seems to 
indicate a crossover between two power-laws (see inset for a logarithmic plot), which occurs around a « 0.6 — 0.7. 
The approximate power-law exponents in the two regions are roughly given by 1.2 and 2.2, respectively. Fig. |5|, on 
the other hand, demonstrates that the dependence of r c on A is almost perfectly linear as expected on theoretical 
grounds and does not exhibit any crossover behavior. The dashed line in Fig. |5| represents the theoretical energy scale 
vTT = ^/ttX. We see that the 20 percent criterion gives roughly r c w 0.7^/Ti. 

The difference between Figs. ^ and || is due to the deviations shown in Fig. || for a > 0.5. This effect is related 
to the finite bandwidth being comparable to the Breit-Wigner width and produces the crossover in Fig. ^. In this 
context, we mention anot her effect which can play a role in such numerical simulations. The mean level spacing of 



the matrix ensemble (1.1) in general depends on a. For large enough values of this parameter the perturbation aHi 
in Eq. (1.1) will start to increase the overall bandwidth of the spectrum and hence the mean level spacing. Therefore 
the dimensionless critical energy r c (a) — u> c (a)/D(a) in general receives contributions from both the numerator and 
denominator. 

In Fig. ^|, we do not observe any crossover behavior showing that r c is a linear function of A. This confirms our 
previous analytical result that the relevant energy scale for deviations from the GUE is given by \/T\ oc A. 
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V. SUMMARY AND CONCLUSIONS 



T he m ain results of the present p aper are the exact expression for A^r) i n ter ms of a handy double integral in 
Eq. fl2.17| ), the perturbative for mula ( 3.31 ) for -X^r) including the extra terms ( 3.42 ), and the perturbative expressions 
for T,' 2 (r) in Eqs. (J0|) and (4.4). Our calculations show that the spectral properties of the Hamiltonian H — Hq + ckHi 
are governed by two energy scales, the Breit-Wigner width T\ and its square root vTT- While the former sets the 
scale for deviations from the Poisson limit, the latter determines the scale at which the difference to GUE statistics 
becomes noticeable. 

To derive these results we employed two supersymmetric approaches, the graded eigenvalue method and the coset 
method. For the technically interested expert, the problem considered in this paper offers the possibility to compare 
both methods in some detail. As a new physical insight from this comparison, we could show that random fluctuations 
of the Breit-Wigner width lead to certain extra contributions to the spectral correlation functions. These contributions 
are automatically included in the (exact) graded eigenvalue method. They are also included in the coset method but 
here the discrete saddle point equation requires a very careful treatment when applying the additional average with 
respect to the diagonal elements because the Breit-Wigner width exhibits small random fluctuations when those 
diagonal elements are varied. 

We briefly discuss the relevance and implications of this paper for previously published work. In Ref. [|l] a dia- 
grammatic analysis of the two-level correlation function X2 (r) was presented. There, a central aim was to p rovide 
a physically more transparent derivation of the results published in Ref. [2^. Comparing our formula ( [2.22 ) with 
the expression (15) in Ref. ^l] it seems that the results of Ref. |2l] are qualitatively but not quantitatively correct. 
This discrepancy, which appears for the second orde r terms, is due to the contributions arising from the statistical 
fluctuations of the Breit-Wigner width (see Sec. Ill E and Appendix ^) which are completely neglected in Ref. |l]. To 
understand this in detail, we have also to take into account that Altland et al. considered a model where the diagonal 
matrix elements of the perturbation aH\ are zero. We have therefore modified their diagrammatic calculation for 
our case w here H\ is a GUE-matrix with non-vanishing diagonal matrix elements. The re sult i s in perfect agreement 
with Eq. (3.31) obtained by the a model perturbation theory. Adding the extra terms (3.41) from the fluctuating 
Breit-Wigner width, we again recover our result (2.22) from the graded eigenvalue method. 

However, we expect that the model of Altland et al. (with H[ 3 = 0) should be equivalent to our case since, in the 
limit DTi/B — > 0, the diagonal elements can be absorbed by an infinitesimal rescaling of Hq- Inspecting Eq. (15) 
of Ref. ^l], we see that here the diffuson perturbation theory already contains the first extra contribution ( |C12|) but 
not the second extra term ( |C25[) . The reason is that for the model with H{ 3 = the fluctuations of Y appear in a 
different way. Indeed, the Dyson (or saddle point) equation is now slightly modified 



T(k) 



1 



DT(j)/2 



(5.1) 



Note that T(fc) depends on the site index and that the sum extends only over the j being different from k. In the 
approximation (|3.13| ) this gives no difference while the trea tmen t of Appendix ^ is modified. It turns out that now 
the first term ( |C12| ) vanishes while the second contribution ( C25 ) still persists. Together with the contributions from 
the diffuson diagrams [Eq. (15) of Ref. 21 1, we obtain again perfect agreement with Eq. ( 2.22| ). 

Apparently, for both models with and without diagonal matrix elements in Hi there are second order contributions 
arising from both the diagrammatic diffuson terms and the fluctuations of the Breit-Wigner width. The sum of 
these contributions is in both cases the same as it should be. However, the way they appear as diffuson terms or 
T- fluctuations is different. 

In Ref. [l^ numerical and analytical evidence was reported for a crossover from a linear to a quadratic dependence 
of the critical energy r c on the strength of thp-perturbing (GOE) matrix elements. From our calculations we cannot 
confirm such a behavior. Both our previousE3 and present analysis indicate the the critical energy r c is linear in 
A, irrespective of the range, A > 1 or A < 1 (see Ref. ^0|), considered. Fig. Ogives the (erroneous) impression of a 
crossover because the strength pa ram eter a is not properly measured on the unfolded scale. The fact that the mean 
level spacing of the Hamiltonian (1.1) depends on a necessitates an a-dependent unfolding. This leads to the linear 
behavior in Fig. 

simulates a very general situation: A diagonalized 



Finally, we would like to re-emphasize that the ensemble (1.1 



quantum system with Poisson statistics is perturbed by a non-diagonal random operator. Situations where this model 
system is appropriate should abound, and in these cases the spectral behavior predicted here, in particular the energy 
scales yTT and Ti should be readily identifiable. 
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APPENDIX A: PERTURB ATIVE EVALUATION OF EQ. (2.19) 



In this Appe ndix we compute the terms of order c 1 and c 2 in the perturbative expansions of the integral P(c, k) 
defined in Eq. (2.19). For this we keep in both sums in Eq. (2.19) the terms up to order p 2 , 



1 /" 27r C°° I P \ ( pe t<fi /pe'^ 2 



-Lq — L\ — L2~ 



(Al) 



with 



2n 



n 



ck 
k + i 



ipr I psmp 



2CK V K 



if n > 1 



(A2) 



Here we have introduced the notion II for the "propagator" of the perturbation theory. To keep a compact notation 
we introduce the abbreviation 



1 



2tt poo 



■(/(./'■ P)){pM = ^J o <W d PP e ~" /(2n) /(^ <P) - 



(A3) 



where /(p, <p) stands for an arbitrary function of p and ip. Expanding the exponential in (Al) with respect to L\ and 
L2, we obtain 



P(c,k) = ( e lip cos 1 



pe -vp , p ,2 



2ni \2ni 



(A4) 



(p,<p) 



Obviously, in (A_4) the terms with odd power in p vanish while the terms with even power can be evaluated using 
{p 2m )(p.ip) = IF"^ 1 2 m ml. The c 1 -term is therefore given by 



P^C^He^cos^W = in 



in agreement with Eq. ( 2.20| ). There are four different second order terms, 



(A5) 



P^M = (^cos^(^) 2 ) = -^cos 2 ^ = - * 



Pf( C , K ) 



e v cos ip — L\ 

ZKl 



1 



4ck 3 



(p 4 cos<^ sin^)( Pj¥> ) = , 



P c (2) (c, k) = -(e lv cos^i 2 ) = -— - (p 4 cos 2 p sin 2 



2!-2 2 n 3 i 



2ck 3 



2ck 3 



it 5 



pf( C ,«)^(e^cos^ 2 ) 
which result in 



1 



' -(p 6 cos 2 ¥ > S in^) = -^-i ? 3!-2 3 n 4 i=- 3 



■ it 



P»fe«>- E ^(,«) = -^(i + 2.(il) +3 (H) 2 



-a.b.c.d 



(A6) 
(A7) 

(A8) 
(A9) 

(A10) 



Inserting II = ck/ (k + i), we obtain Eq. (2.22) 



1G 



APPENDIX B: PERTUBATION THEORY IN THE FRAME WORK OF THE COSET METHOD 



The precise a model action Cq and the preexponcntial terms in (3.21) still depend on the particular realisation of 



the diagonal matrix elements Hq J . The aim of this appendix is to show how one can perform in a suitable way the 
average with respect to these matrix elements and how the resulting Q-integral can be evaluated perturbatively for 



Q being close the origin A in the coset space. In order to expand Cq with respect to AQ 
consider the quantity 

/ = trg ln(a + bA + cQ) 



{A,Q}/2- 1 = 2R 2 , we 



(Bl) 



with some complex parameters a, b, c and A, Q as defined in section III C . The quantity / only depends on the radial 
parametersa 8 of Q — T _1 AT, i.e. n we may choose T = exp(0<Ti), where o\ = L J) and 9 is 2 x 2 supermatrix with 
two independent radial parametersu (for the unitary symmetry class). Then we have a\Aa i — —A and <T\Q<J\ = — Q. 
Inserting 1 = u\ and permuting the matrix products in the argument of the logarithm in (|Bl|), we thus obtain: 



I = trg ln(a - &A - cQ) = -trg ln( (a + bA + cQ)(a -bA- cQ) 



—trg In 
2 B 



a 2 -(b 



-be 



= 2 tr S ln ( 1 



26c 



(b- 



QA 

aq" 



Using this identy, we can rewrite £ as 



AQ-2 



= —trg In 

2 6 



{b + cf -2bcAQ 



(B2) 



4 = ^trg In 1 



iuoTiD/2 



Htf + T D/2) 2 + (Tj.D/2 - ioj/2)< 



AQ 



(B3) 



(2) t 

The logarithm can be exanded in powers of AQ giving £ = £ y Q ' + C, 



,(4) 



where Cq is the gaussian part of the 



action. Furthermore we note that in the limit Ti 3> 1 the discrete j-sum consists of Lorentzians with well overlapping 

— 1 (2) 

resonances which allows to replace in lowest order of r 1 the gaussian part Cq by its average over the diagonal 
Therefore we obtain, in the limit B/Y\ — > oo, 



elements H J 3 



(2) 



- % n r n ■ tr S AQ =' l o r — T" trg AQ 
D L iD — iuj 2 1 1 — ir 



(B4) 



where r — lo/D and B is the bandwidth over which the diagonal matrix elements Hq 3 are distributed. Omitting the 
details, we mention that one can show in a more refined way that the statistical fluctuations indeed yield corrections 

— 1 (2) 

of lower order in T x . For this, one has to analyze the average of the quantity exp(— Cq (Q)) more carefully. 

The spirit of the pertubative calculation consists of expanding the the preexponcntial terms [and also of the non- 

gaussian terms exp(— — • • •)] with respect to the variables R (or t) [see ( [3.25 )]. Then one can systematically 



evaluate the gaussian integrals using Wick's theorem. Each pair of t and t yields here a power II^, 1 of the propagator 



which is defined by the prefactor of the gaussian action 



4\Q) 



1 1 ' trg(i? 2 ) = -2ni J Tl . trg(tt) = U D trg(tt) 



Ti — ir 



Ti — ir 



(B5) 



The perturbation theory is justified for |IId| ^> 1 which is certainly correct if Ti ^> 1 and if the ratio k — r/Y\ is 
finite and non 0. As in the pure GUE case, this simple perturbation theory does not allow to access th e ver y small 
frequencies to ~ D. We emphasize that this is different for the pertubative evaluation of the integral ( 2. IE ) which 
was obtained in the frame work of the precise graded eigenvalue method in section [H] (see appendix |A). Here the 
pertubative cond ition only consists of T\ 3> 1 and the behavior for small u) is precisely described by the residuum 
contributions in (2.16) which even contain the nonanalyti cal co rrections. 

We have to expand gj in the preexponential factors in (3.21) with respect to the small quantity 



Q - A = 2Ai? v / l + R 2 + 2AR 2 = 2AR + 2AR 2 



(B6) 



which gives 
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(0) 

9 3 = 9) 



r i-D\ (o)^ (o) 
«— — I g) SQg) 



where 



(B7) 



(0) 



,(0) 



x - H 3 j + T D/2 + iA(T 1 D-iuj)/2 
x - H 3 j + T D/2-iA (Ti D - iuo) /2 



To evaluate the term Si(Q) in ( |3.2lD , we consider 



^trg( ft P ± )=^trg( 5 f P ± ) 



riD 



]>>g[(. 9 f ) 2 (2AP + 2AP 2 )P±] 



(B8) 
(B9) 

(BIO) 



+ Y,^\ 0) ) 2 gf ) (2AR2AR)P ± } 



.(±) 



.4 



(±) 



.4 



(±) 



where we have used that g^ R — R g^' since APA = — R. The term Si(Q) is now given by 



5i(Q) = + A +) + Af>) I A''' + A\-> + A [ 2 



i + A (aV 



(Bll) 



The main contribution arises from the product A^ Aq ^ which is effectively of order II^, and also contains terms of 



higher orders IT^ 1 , rip 2 , .... This contribution just corresponds to the boundary terms ( 3.22 ) which are extensively 
discussed in subsection HIE and appendix |c[ Here we concentrate on the higher order terms arising from A\ 
and A2 which are both of second order in R (the linear R term in A\ vanishes under the graded trace). They 
will give contributions of order II^, 2 which is the highest order in which we are interested. (This corresponds to 

order c 2 in appendix 1X1.) Therefore we can here replace the correlated Po-averages like [A2A2 exp(— Cq )]q by 
(2) 

[A2]o [A2]o exp(— [£q ]o) since the correlations arising from the diagonal terms or the exponential are already of order 
IT^) . The crossed terms like [j4o^2]o ar e in principle of order II^ 1 but they vanish upon the gaussian contraction 
rules (for details see Ref. ||^). Actually, they are of the type trg(iiP+) — > II^ 1 trg(l) trg(P+) = 0. The product of 

these terms with first term arising form the nongaussian part of the action, i.e. [Aq A2\q£ < ^\ is of order II^ 2 and of 
the type trg(itP+) tig(titt) which also vanishes. The term [^4i]o [Ai]o vanishs upon the Pp-average while [A 2 ]o [^b]o 
gives the only nonvanishing contribution. The Po-average can be performed using the integrals (C6,C^) evaluated in 
appendix |c]: 



A 



(±) 



TxDy {-2m) 



D {TiD-iuj)' 



trg(-4AP 2 P±) 



2?ri 



Ti — ir 



trg(i iP+) 
-trg(t*P_) 



In summary, we obtain for the term Si(Q) 

[5 1 (Q)]o»5 1 (A) + ^ 



tvg(ttP+)trg(ttP-) 



(B12) 



(B13) 



which is just the first part of Eq. (3.27) 



Concerning the term ^(Q), only the linear SQ-terui in (B7) contributes 



S 2 (Q) 



^ trg (2AP) gf P + gf (2AR) gf P_] 



(B14) 



£trg[( 5 f) 2 (gf) 2 (-4)PP + PP_] 



We can again perform the Pp-average by using ( |C7|]C8| ) 
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[S 2 (Q)h 



1 



HO 7; 



4tt 



D (TiD-iu)- 



4tt 2 

trg [PP+PP_] = -j. 



P> 2 (r!-ir) ; 



■trg[f P+tP_] 



This expression provides the seco nd p art of Eq. (3.27). 

The gaussian contractions for (B13) and ( B15| ) can be performed by 



trg(«P+)trg(UP_) n^trg[i P+tP_] n^ 2 trg(P + ) trg(P_; 



n 



and 



trg[t P+tP-] -> n^trg(P + )trg(P_) = n 



(B15) 



(B16) 



(B17) 



since trg(P±) = 1. Combining this with (B13) and (B15), we finally obtain the result ( 3.28| ). 

We mention that this perturbative calculation can be immediately generalized to the orthogonal (/3 = 1) and 
symplectic symmetry classes ((3 = 4). In these cases t is a 4 x 4-matrix containing 2 x 2-blocks and tg for the 
diffuson and cooperon contributions. The result of this generalization just gives the prefactor 2/(3 in ( 3.28| ). 



APPENDIX C: CONTRIBUTIONS FROM THE FLUCTUATING BREIT— WIGNER WIDTH 



The statistical fluctuations of the Breit-Wigner width V discussed in subsection [HE result in further second 
order contributions when the boundary term S±(A) is averaged over the diagonal matrix elements. To extract these 
contributions, we need the following type of simple integrals 

Wt) = L dy (^TW (v - *W • (C1) 

Obviously we have 

/i,0 = = -«r , I n ,o = lo.n = if n > 2 , (C2) 

because for n > 2 the integration contour can be closed such that there are no singularities inside the contour. For 
77i = 1, the pol of first order inside the contour results in 

/„,i(7) = H)"- 1 % , n>l. (C3) 

Furthermore, we obtain by partial integration the recurrence relation 

71 

In,m{l) = r I n +l, m -lh) (C4) 

777 —1 

which can be solved in terms of the binomial coefficient 

InMl) = (-l) n - 1 i n+m - 2 ( n l l T 1 2 )y^ > ^,m>l. (06) 
For the subsequent calculation, we need the cases 

h,i{-y) = — , hail) = -hAi) = i% , (C6) 

7 7 

/i,a(7) = ^3,i(7) = , ^,2 (7) = 2 % . (C7) 

ryJ ry<J 



The moments of the quantities Gj, Gj [see Eq. ( p.3q )1 with respect to the H average can be expressed in terms of 
these integrals, 



_ 1 i-b/2 1 1 1 

2J G " G f] = d b^Ic J_ B/2 dH [(x-H)+ 7(^-7^/2]" [(x -H)- i{DTi - iu)/2\™ = D In ^ DTl ~ ^ 



(08) 
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Now, we turn to the average of the boundary term (3.4C) which we rewrite as 

Si (A) = (B + B 1 )(B + Bi) 



(C9) 



with B = Y^,j Gj, B\ = —D5T/2J2jG^ and analog for B , B\. The main contribution [B B } has already been 
considered and evaluated in (3.23). The next contribution is given by [see (3.37)] 



j,k j,k 



Gj Gfc 



D 



E G 



(CIO) 



DTt 



E 



G ? g j ( ~ E G 



The terms with k ^ j vanish since all averages of the type [G™]o, [G™]o vanish if n > 2 [see ( |C2[ )]. Furthermore, 
the contributions with I ^ j cancel for N — > oo the second term because X)/[Gz]o = Ii,o(DT\ — iu>)/D = —iir/D. 
Therefore, by (^|@ we obtain 



[**]o = ^£[<^] Z)2ri 1 



2tt 



2tt D V (-CTi - iw) 3 / L> 2 (fi - ir) 



i ri 



In an analog way, one finds that [Bq B±]q gives the same contribution, such that 

2 r x 



[B\ B ] + [B Q Bi] 



The remaining contribution, 
D 2 T 2 

[BiSi]o = ^i 



E 



✓~<2 /^i2 



L> 2 (Tj - ir f 



E G <+0(?E G 



C4 + C5 + C6 



leads to three different subcases with 4, 5, or 6 G-factors. In the first case only the terms j = k contribute 

n 2 r 2 ^ D 2 r 2 



Ci - E [ g j G ^] ~ E [ g j G T\ 



(Cll) 



(C12) 



(C13) 



(C14) 



For G5 and G«, we have to consider all different cases of equal or different summations indicies. In the following, we 
mention only those cases giving a nonvanishing contribution. Actually, for 



' E [G 2 J Gt(G m -G m )] Q 



4tt 



there are two relevant cases: 

G5 : j = k^m 

G 5 : j = k = m 



n 2 r 2 

-^E^o 



3^2 



d 6 t 

47T 



i E[ G l G K G J '- G i)], 



The last term 



Gk = 



.p 4 r 2 

4tt 2 



E [ G ? G fe Gl Gm \ 



(C15) 

(C16) 
(C17) 

(C18) 



leads to six cases: 
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C 6 : l^(j = k)^m^l 

C 6 : (j = k)^(m = l) 

C 6 : (j = m)^(k = l) 

C 6 : (j = k = m)^l 

C 6 : (j = k = l)^m - 

Cq : j = k = I = m — > 



D 2 T 2 



E [ G i G T\ o 



4tt 2 

D 4 T 2 
4tt 2 



X! t G ? G ?] o X! [ Gm Gm ) o 

j m 

H [° 2 j G J'] E [ Gfe G fc] 



i>T 2 

47T 



E [ G i G i] o 



3p2 



z>T 



4tt 



z? 4 r 2 - i r 2 
-± J] [g 3 - g 3 ] - — _\ - (n z 



4tt 2 ^ " JJ Q " £> (Ti -ir) 5 ^ 1 ' 7J 



(C19) 
(C20) 
(C21) 
(C22) 
(C23) 
(C24) 



We observe that (|C14|), QClg ) and ( )C19| ) mutually cancel. Furthermore, flC17| ) is cancelled by (C22J) and (|C23|) ( |C24| ) 
is already of third order and has to be disregarded. The only second order contributions are ( p20| ) and (C21) which 
give by (p3-[C8|) 



[B! Bi] 



£> 4 F 2 1 



4tt 2 



o r 2 

[2 + i(-i)] = -^ 



4tt 2 D 2 (DTi - iuj) A L " ' " v D 2 (T 1 - ir) 4 ' 

This expression and (C12j) are just the extra terms mentioned in (3.41). 
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FIG. 1. Graphical illustration of the shifts of the complex integration paths. The two integration paths (full lines) differ by 
the closed contour integral (dashed line). The direction of the shifts are indicated by the two thin-dashed lines. The first shift 
2/i = xi +i i s show n in (a). Here for \y2\ < 1 the pole — i + |j/2| lies inside the contour and results in the residuum contribution 
_Ri [see Eq. (2.17)]. The second shift 112 = y2 + in is illustrated in (b) (for k > 0) and (c) (for k < 0). T he singularity ±i|i/i| 
lies inside the contour if < \n\. The corresponding residuum contribution is given by R2 in Eq. (2.17). 
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FIG. 3. The parameter A as a function of a. Shown is the value of A obtained from the fit of our analytical formula for the 
number variance to the simulated data (full line and points) and the theoretical expression A th = a/D (dashed lined). Here D 
is the mean level spacing for the parameter values mentioned in the text. Note that the critical value a c ~ 0.95 corresponds to 
the situation in which the Breit-Wigner width Fi D equals the radius Ro of the unperturbed semicircle density. 
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FIG. 4. The dependence of the critical energy r c as defined by the 20 percent criterion explained in the text on the strength 
parameter a. Inset: The same on a logarithmic scale. In addition we have drawn two dashed lines with slopes 1.2 and 2.2, 
respectively, to guide the eye. 
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